Full frequency voltage noise spectral density of a single electron transistor. 
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We calculate the full frequency spectral density of voltage fluctuations in a Single Electron Tran- 
sistor (SET), used as an electrometer biased above the Coulomb threshold so that the current 
through the SET is carried by sequential tunnel events. We consider both a normal state SET and 
a superconducting SET. The whole spectrum from low frequency telegraph noise to quantum noise 
at frequencies comparable to the SET charging energy (Ec/h), and high frequency Nyquist noise is 
described. We take the energy exchange between the SET and the measured system into account 
using a real-time diagrammatic Keldysh technique. The voltage fluctuations determine the back- 
action of the SET onto the measured system and we specifically discuss the case of superconducting 
charge qubit read-out and measuring the so-called Coulomb staircase of a single Cooper pair box. 

PACS numbers: 03.67.Lx 42.50.Lc 73.23. Hk 85.25.Na 



I. INTRODUCTION 

Solid-state realizations of qubits are of real interest due 
to the possibility of using lithographic techniques to in- 
tegrate the large number of qubits, needed for a fully 
functional QC. The Single Electron Transistor (SET) has 
been suggested as r%-|Hftfl" ou t device for different solid- 
state charge qubitsCJBBcl'B. 

Aassime et al.H have shown that the Radio-Frequency 
SETS (RF-SET) may be used for single shpt read-out of 
the Single Cooper-pair Box (SCB) qubitoH. This is pos- 
sible if the measurement time t ms needed to resolve the 
two states of the qubit is much shorter than the time t m i x 
on which the qubit approaches its new steady-state deter- 
mined by the back-action due to voltage fluctuations on 
the SET. In a previous short paperfj we provided further 
support for this result by calculating the full frequency 
voltage noise spectral density of the SET, including the 
effect of energy exchange between the qubit and the SET. 

In this paper we give a full account of the calculation as 
well as a thorough discussion of the effect of back-action 
in measuring the so-called Coulomb staircase of an SCB 
qubit. We also include a section about the back-action 
from a superconducting SET. 

Since the qubit is carefully shielded from all unwanted 
interactions with its environment it is reasonable to as- 
sume that the back-action from the SET charge measure- 
ment is the dominating noise source, even though the 
two systems are only weakly coupled. This further moti- 
vates choosing measurements on a charge qubit to discuss 
the scectral properties of SET back-action, compare e.g. 
Rcf.EI. In non-qubit systems, e.g. a normal state single- 
electron box, other sources of dissipation dominate over 
the back-action from the charge measurement. Further- 
more we do not discuss the deph asing of the qubit in- 
duced by the presence of the SETliilii3. Although it is a 
very important subject, the dephasing time is mainly de- 
termined by the zero frequency fluctuations, and in this 



limit our result coincides with previous expressions!!!! 

The structure of the paper is the following: In sec- 
tion H we discuss the basic properties of the SCB qubit 
and the effect of gate voltage fluctuations. Furthermore 
the difference between asymmetric and symmetric def- 
initions of noise spectral density is noted and also the 
connection to mixing-time and the Coulomb staircase. 
In section III the model for the SET is introduced and in 
section IV the real-time diagrammatic Keldysh technique 



is described and the expression for the spectral density of 
voltage fluctuations on the SET island is derived. Section 
^ describes the properties of the voltage fluctuations in 
different frequency regimes, the effect on the mixing-time 
and the Coulomb staircase both for a normal state SET 
and a superconducting SET. 



II. THE SINGLE COOPER-PAIR BOX QUBIT 

The qubit is here made up of the two lowest lying en- 
ergy levels in a single Cooper-pair box (SCB)I3. An SCB 
is a small superconducting island coupled to a supercon- 
ducting reservoir via a Josephson junction. The Hamil- 
tcuaian of the system can be written in the charge basis 

bJb 

H q = J2^E qb (n-n g ) 2 \n)(n\-^- [\n)(n + l\ + \n+l)(n\] 

n 

where E q b = e 2 /C q i> is the charging energy, n is the num- 
ber of extra Cooper-pairs on the island, Ej is the Joseph- 
son energy of the junction and n g — Cf'V^ /2e is the 
number of gate-induced Cooper-pairs. In order to get a 
good Cooper-pair box we need ksT <C Ej <C E c < A, 
where A is the superconducting gap and T is the tem- 
perature. The low temperature is required to prevent 
thermal excitations and the high superconducting gap is 
needed to suppress quasiparticle tunneling. The eigenen- 
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FIG. 1: The energy bands of an SCB. The states |0) and 
|1) denote the eigenstates on the island. The quasiparticle 
branches have been left out for simplicity. The gap is due to 
hybridization of the charge states by the Josephson coupling. 



ergies now form parabolas, varying with the gate voltage. 
For suitable values of the gate voltage (close to n g = 1/2), 
and for Ej <C E c , the system reduces to an effective two 
level system as the two lowest lying charge states are 
well separated from the states with higher energy. In- 
cluding only the two lowest lying states in Eq. (]]]) the 
qubit Hamiltonian becomes 

H q = -^(l-2n g )a z -^-a x , (2) 

where a x . z are the Pauli matrices (and the states |j) = 
{ P and \l) = ^ ^ ^ correspond to zero and one ex- 
tra Cooper-pair on the qubit island). By changing the 
gate voltage the eigenstates of the qubit can be tuned 
from being almost pure charge states to a superposition 
of charge states. 

The ground state/first excited state of the system, 
written in the charge basis [Ej — 0) are 

|0) = cos(V2)|T)+sin(7 ? /2)|I), 

|1) = - Bin(tj/2)|T) + cos(r ? /2)||), (3) 

where r\ = arctan(£ , j/4£' 9 ( ) (l — 2n g )) is the mixing an- 
gle. The energy difference between the two states is 
AE = ^(4:E q b) 2 (l — 2n g ) 2 + Ej and the average charge 
of the eigenstates is 

Q =2e|a |0)| 2 = 2esin 2 (r y /2), 
Q 1 = 2e|(j |l)| 2 = 2ecos 2 (r7/2). 

(4) 



A. Qubit transitions induced by SET voltage 
fluctuations 

When the SET is turned on in order to measure the 
qubit, the voltage fluctuations on the SET-island induce 
a fluctuating charge on the qubit island. This is equiva- 
lent to a fluctuating qubit gate charge n g — > n g + 5n g (t) , 
giving rise to a fluctuating term in the qubit Hamilto- 
nian, written in the charge basis as 

AE h 

SH q (t) = —^L25n g {t)a z = 2e K SV(t)a z , (5) 

where k = C c /C q b- Here SV(t) represents the voltage 
fluctuations on the SET-island, and we have neglected a 
term quadratic in 5n g (t). In the qubit eigenbasis, using 
the rotation defined by Eq. (||) , the fluctuations in Eq. (g) 
become 

5H(t) = 2en8V(t) [cos{rf)a z + sm(r))a x ] . (6) 

The fluctuating voltage on the SET island can induce 
transitions between the eigenstates of the qubit. If the 
capacitive coupling to the SET is small (k <C 1) we can 
use the Fermi golden rule to calculate the transition rates: 

T rel (AE) = *J^ K *S v (*E/h), (7) 

T exc (AE) = ^-^LtfSvi-AE/h), (8) 

where T re i is the relaxation rate and T exc is excitation 
rate and Sy(AE) is the asymmetric (see Eq. ([)])) spec- 
tral density of the voltage fluctuations on the SET island. 
The fraction Ej / AE — sin(?7) comes from Eq. (j|) as it is 
only a x that causes any transitions between the states. 
Note that the cos(rj)a z term causes fluctuation of the en- 
ergy levels, leading to phase fluctuations and dephasing. 

B. Asymmetric noise - Coulomb staircase 

In our calculations we emphasize the energy exchange 
between the qubit and the SET and separate between the 
contributions from processes leading to the qubit loosing 
energy and the contribution from processes leading to the 
qubit gaining energy (or equivalently: the SET absorb- 
ing or emitting energy). Because of this, we will maintain 
this separation of the noise spectral density of the SET 
into contributions from positive and negative frequencies 
and therefore use the asymmetric expression for the volt- 
age fluctuations 

/oo 
e~^ T (5V(T)6V(0)). (9) 
-oo 

As we are primarily interested in the processes in the 
SET we chose our reference so that positive frequencies 
correspond to the SET absorbing energy and negative 
frequencies correspond to the SET emitting energy. 



One example where the separation is necessary is in 
describing the back-action of the SET while measuring 
the so-called Coulomb staircase, i.e. the average charge 
of the qubit as a function of its gate voltage. In an 
ideal situation with no energy available from an exter- 
nal source, at zero temperature, the qubit would follow 
the ground state adiabatically and the charge would in- 
crement in steps of 2e at n g = n + 0.5, n integer. These 
steps are not perfectly sharp because of the Josephson 
energy mixing the charge states. This mixing of charge 
states results in a maximal derivative given by 4Ec / EjU. 
When adding a noise source, e.g. by increasing the tem- 
perature or attaching a noisy measurement device, the 
steps will be rounded further due to a finite population 
of the excited state. 

Assuming that the SET is the dominant noise source, 
and that the two state approximation of the qubit is valid, 
the steady-state population of the qubit is given by 




FIG. 2: The SET. 



III. SET MODEL 



We follow the outline of Ref and model the SET by 
the Hamiltonian 



^st^qb 



T exc (AE)/[T exc (AE) + T re i(AE)} (10) 



P 



st,qb 



T rel (AE)/[r exc (AE) + T rel (AE)}. (11) 



H = H L +H R + H I + V + H T = H + H T , (15) 
where 



Hr 



e kn a krn a km, #J = E e '« c L c i™ ( 16 ) 



The corresponding expression for the average charge is 
then given by 

Q(AE) = Q q (AE)P^\AE) + Q 1 (AE)Pf' qb (AE), 

(12) 

where n is the charge of the excited state/ground state 
defined in Eq. (§). If we set Pf' qb = and P° t ' qb = 1, we 
recover the ideal Coulomb staircase, as this corresponds 
to the system following the ground state adiabatically. 



C. Symmetric noise - Mixing time 

For quantities that depend on the summed rate of 
relaxation- and excitation-processes in the SET the sepa- 
ration of absorption and emission might not be necessary, 
and the symmetrised expression for voltage fluctuations 



describe noninteracting electrons in the left /right lead 
(H r ,r £ {L,R}) and on the island (Hi). The quantum 
numbers n denote transverse channels including spin, and 
k, I denote momenta. The Coulomb interaction on the 
island is described by 



(17) 



V(N) = E c (N-n x ) 2 



where N denotes the excess number op- 
erator, Ec = e 2 /2C the charging energy 
(C = Cl + Cr + C g + C c ), n x the fractional num- 
ber of electrons induced by the external voltages (n x is 
the fractional part of (ClVl + CrVr + C g V g )/e) and e 
the electron charge. The tunneling term is 



Hi 



E 



L.R kin 



,ci n e 



rprn* T 
1 kl c h 



O-kr 



H~ +H+, 



) 

(18) 



Sy ym (u) = / dTe-^ T (SV(T)SV{0) + SV(0)SV(t)) 

J — oo 

(13) 

can be used. Note that Sy m (uj) = Sv(w) + Sv{—uj). 

One example is the time it takes the qubit to reach 
its steady-state, after the SET is switched on. This time 
is called the mixing time, and the information about the 
initial state-nonulation is lost on this timescale. For weak 
coupling it isLftLS 



r reZ (A£) + r e 



,(AE) oc Sy ym (AE/h). (14) 



where the operator e ±J * changes the excess particle 
number on the island by ±1 and T£™ are the tunnel- 
ing matrix elements. <£ is the canonical conjugate to 
N, ([<t,N] = i). In this case of a metallic island con- 
taining a large number of electrons, the charge degree of 
freedom N — 0, ±1, ... is to a very good approximation 
independent of the electron degrees of freedom I, n. The 
terms H + and H~ represent electron tunneling to and 
from the SET island. The form of the tunneling terms 
(with the e** term) is a consequence of separating state 
space into electron /, n and charge N degrees of freedom. 
This is also reflected in the partitioning of the density 
matrix introduced below. 
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The spectral density of voltage fluctuations on the 
SET island is described by the Fourier transform of the 
voltage- voltage correlation function 



the equilibrium (Fermi distributed) density matrix p e eq 



for 



dre 



e 

C 2 



= t^t / dre 



r Tr 



r Tr 



{p st (t )SN{T)6N{0)} = 

[p st (t ) (iV(r)7V(0)-iV 2 )} 

(19) 



the electron degrees of freedom in each reservoir (L, R, I) 
and a reduced density matrix p c st , describing the charge 
degrees of freedom. Since the tunneling events between 
the SET island and the electrodes are incoherent due to 
the low conductance of the tunnel junctions, the charge 
states wilLbe incoherent, and p c st is therefore taken to be 
diagonals with elements Pj$ denoting the steady-state 
probability of being in charge state N. 



where TV is the average of N. The /V 2 -term in Eq. 
assures that the correlation function vanishes for large r. 
In Fourier space this term does not contribute at finite 
frequency, and at zero frequency it compensates for the 
steady-state delta function. In order to simplify our ex- 
pressions we leave this term out, and keep the frequency 
finite during the calculations. 

In Eq. ( JlOj ) p s t(to) is the density matrix of the sys- 
tem in the steady-state, which is assumed to have been 
reached at some time to before the fluctuation occurs 
(to < min{0, r}). p st = p e eq ® p c st is the tensor product of 



IV. DIAGRAMMATIC TREATMENT 

We now expand the correlation function in Eq. ( |l9| ) in 
a perturbation series in terms of the tunneling Hamilto- 
nian (Ht) along a Keldysh time contour (see Fig. ||). The 
trace over the electron degrees of freedom is then eval- 
uated using Wick's theorem, which is possible since the 
tunneling Hamiltonian is only bilinear in the fermionic 
operators. Rewriting Eq. in the interaction picture 
gives (excluding the iV 2 -term) 



SV(u>) 



C 2 



dre 



jr Tr {p(t )S(t , t)N(t)S(t, 0)N(0)S(0, t )} , 



where S(t%, ti) is the S'-matrix that brings the system from the time t\ to time ti^ i.e. for t2 > t\ 



,M) = e dtHT(t) = l-i T dnffrCn) + H) 2 r dn r dr\ dr2 Ht (j\ ) Ht ( t 2 ) 



(20) 



(21) 



and analogously for i 2 <t\. Defining T as the time-ordering operator along the Keldysh contour we can write Eq. (EC 



„2 poo 



dre 



Tr 



{p(t o )N(r)N(0)} 



driTr 



{T jO (t )7V(r)iV(0)ffT(Ti)} 



+ 



H) 

2! 



2 r ° r° r i 

dr x / dT^\Tp(to)N(T)N{$)HT(T X )HT(T2)\ 



+ . 



(22) 



r 



Since the unperturbed Hamiltonian does not contain 
any couplings between the leads and the island nor be- 
tween the leads themselves, their degrees of freedom are 
independent. Moreover, as the trace of independent de- 
grees of freedom is equal to the product of the respective 
traces and every perturbation term contains one operator 
from one of the leads and one operator from the island, 
only terms containing an even number of perturbation 
terms will contribute. 

Using the diagrammatic language of Ref.0, the noise 
correlation function SV(w) can be given a diagrammatic 
formulation. 

Fig. ^| shows a diagrammatic representation of the first 
term in Eq. (E3) with zero Ht perturbations, involving 



the evolution of the state described by the density ma- 
trix itself along the Keldysh time contour. This process 
is represented by a solid line starting and ending at the 
the density matrix, including the charge (fluctuation) op- 
erators N(t) and iV(0) marked by dots. This gives the 
expectation value (statistical average) of the charge fluc- 
tuation correlation (in the " ground state" ) in the absence 
of any transport process. 

The first contributions to charge transport come from 
the third term in Eq. (|22|) with one H + and one H~ per- 
turbation, describing tunneling onto the SET island and 
back, changing the charge state N to N ± 1. A diagram- 
matic representation of this would be the middle part of 
Fig. ^. Since the charge-transfer process can be viewed 
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N(0) 



N(T) 

— I 1— 1 



FIG. 3: The propagation of the unperturbed steady state 
density matrix. This corresponds to the first term in Eq. 



N(T) 




1 1 — t 

T 

FIG. 4: An example of a specific diagram. The wiggly lines 
correspond to island lines and the solid lines correspond to 
lead lines. The rightmost part correspond to an electron-hole 
excitation while the leftmost part corresponds to a correction 
to the external vertex. The middle part corresponds to a 
tunneling event as only the diagrams connecting the upper 
and lower branch changes the number of extra charges on the 
island. 



as an electron-hole excitation, creating a hole on an elec- 
trode and an electron on the island, in Fig. |] there are 
new additional lines with arrows representing electron- 
hole propagators (excitations). Every internal time will 



form a vertex and the propagator (Ta(Ti)a 1 ' (t 2 )) will 
form a line going from r 2 to T\. In this case with macro- 
scopic metallic reservoirs withjnany transverse channels, 
the main contributing termalj will appear in combina- 
tions of 

(th+( T1 )h-(t 2 )) = EE E E [nirnT * 

ri r 2 kxhni k 2 hri2 
X<^a*irini(Tl)4 a r 2 n s ( r 2)) X 
X ( T 4 iri n 1 ( T l) C fer 2 n 2 ( 1 '2)) , 

(23) 

which means that there will always be pairs of internal 
lines with reversed start and end points, one being a 
reservoir propagator and one being an island propaga- 
tor, as shown in Fig. 0. These line-pairs are replaced 
by a single line, corresponding to an electron-hole exci- 
tation, with an energy equal to the difference in energy 
between the two lines. 

In order to facilitate the evaluation of the time-ordered 
diagrams in Eq. (122) , we rewrite the Fourier-transform in 

e = e 2 -e 1 

FIG. 5: Two internal electron lines with reversed start and 
end points. They can be replaced by a single line with an 
energy equal to the difference in energy between the two cor- 
responding to an electron-hole excitation. 



Eq. (||) as 



S v (uj) 



dre-^Tr{p at (t )N(r)N(0)} = — 



' — oo 
„2 r 



dTe- iuT Tr{p st (t o )N(r)N(0)} 



C 2 



/ dre+^ Tr{p st (t )N(0)N( T )} = 



dre- iUT Tr{ Pst (t o )N(T)N(0)} 



(24) 



where we have used that the steady-state is time invari- 
ant. Furthermore we fix the specific time ordering of all 
internal and external times in all diagrams. This make 
the diagrams straightforward to evaluate in the frequency 
domain as all integrals thus become recursive Laplace- 
transforms. 



Returning to Eq. (^2|) and drawing all the diagrams of 
the lowest non-trivial order we can divide them into two 
categories: Dressings of the propagator H N ^ N >(cu) that 
takes the system from the charge state N to the state N' 
(see Fig. ^), and vertex correction (see Fig. 0). 



Calculating for instance the diagram in Frig. ^ using where An — V(N + 1) — V(N) is the charging energy 
the diagrammatic technique outlined in Ref.li-3 yields cost of moving one electron from one of the leads onto 

7+(£) 



D N = -Y lim - / dE- , (25) 
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N+l N-l N-l 




N N N 




FIG. 7: All the vertex corrections of the external vertex N(t), 
denoted by a dot, within the sequential tunneling approxima- 
tion. Note the correspondence to Fig. pi 



the island when there are N extra electrons there. The 
factor 7^(e) is the inelastic (Golden Rule) tunneling rate 
through junction r for electrons tunneling to the island 
(correspondingly 7" (e) is the inelastic tunneling rate of 
electrons tunneling from the island through junction r) 
given by the expression 



^( £ ) = * E / dE P( £ + E ~ eV r )p(E)f ± (e + E- eV r )F{E)\T™\\ (26) 

I 



where p is the density of states in the reservoir r and 
on the island, both assumed to be either superconduc- 
tors or normal metal, f + (E) is the Fermi-distribution, 
f-(E) = 1- f+(E), T rn is the tunneling matrix ele- 



ment, here assumed to be independent of energy and 
V r = (p r — Pi)/c- is the voltage bias across the r junction. 
Separating Eq. ( p5| ) into real and imaginary parts we get 



lim - / dE~f+(E) 



Tiuj + E - A 



II] 



(huj + E- A n ) 2 +t] 2 (hu + E-A 



(27) 



The real part of the integral is small in the tunnel limit malization effects (see Ref.E-3) and concentrate on the 
where the conductance is small, so we neglect these renor- 
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imaginary part, which gives 

r 

Introducing the notation 



r=B.,L 



r=R,L 



«)^[7^(0)-7^M], 



(28) 



(29) 



d) ^[7^(0) 

e) ^(0) 



all the diagrams in Fig. |6| can be calculated in the same 
way, resulting in 



a )^"7wM: c0^7atM, ff)=*"-7jvH 
6)^ 7 +(-o;), e)=>7w(-w), M =* (30) 
c) -7jv(-^), /) -7iv(-w) i) 

In the same way, the vertex corrections in Fig. can 
be calculated yielding 



7tv(-^ 



<0 - ^te(O) - 7* ("")], /) - i[7j^(0) - 7^(- 



9)-^[7^(0)-7^M] 
> h)-±tr+(p)-<Y+(w)] 
oj)}, 



(31) 



where the expression in Eq. (|31|i) corresponds to the ze- 
roth order correction of the vertex. To lowest order, the 
total spectral density can be written 



S v (u)) 



2 h Re 



E 



N,N',N" 



2^ Re \^ T fl(cu)V(iu)V St } 



(32) 



where ~$ and ~P are column vectors is the trans- 
pose of JV) containing the number of extra electrons on 
the island and the steady-state probabilities respectively. 
Note thus that the element Pq 1 refers to the steady-state 
probability of the SET to be in the state with ex- 
tra charges, unlike for the qubit, where p^' qh refers to 
the steady-state probability of being in the lower energy 
eigenstate. V(u>) is a matrix whose elements Vn>,n(w) 
are the sum of all vertex corrections in Fig. [?] which take 
the system from the state N to the state N' and II(w) 
is a matrix whose elements ILn",n'(oj) are the sum of all 
diagrams in Fig. ^| that take the system from the state 
N' to the state N" . Incidentally, Eq. @ is valid for 
arbitrary order, as long as the vertex correction and the 
propagator are dressed to the appropriate order. 

In order to facilitate the evaluation of higher order di- 
agrams we define an irreducible diagram as a diagram 
where it is impossible to draw an auxiliary vertical line 
at any time, without crossing an electron-hole line. An 
example can be seen in Fig. |S|. 

The diagrams in Fig. Hare all the first order irreducible 
diagrams, except Fig. |6j), which is just the free propa- 
gator flff N ,(uj) = ^5n,n', where Sn.n' is a Kronecker 
delta. 

Using irreducible diagrams allows us to write down a 
matrix Dyson equation in frequency space for the fre- 
quency dependent propagator II (u>) between different 
charge states (Fig. 0) 



ZD 



FIG. 8: This diagram contains three irreducible diagrams 
as in between it is possible to draw a vertical auxiliary line 
(dashed lines) separating them. 



N' N 



n«A) 



n>) 



N' N 



+ £ 

.. N",N" 



lT(ffl) 




E(w) 




new) 


N,N" 




N* %N'" 




N"',N' 


N 


N 





FIG. 9: A graphical representation of the Dyson equation in 
Eq. (Q). The terms T, NN i in are irreducible diagrams that 
take the system from the state N to the state N' . 



ft(w) = n (0) H + nW(w)E(w)n(w). 



(33) 



Solving for II (u) and inserting the explicit form of 
n(°) (uj) we get 



1%) 



(34) 



Note that the explicit time ordering in every diagram, 
means that in frequency space any specific diagram can 
be written as a product of irreducible diagrams and free 
propagators. 

The matrix element Sjv',j\r(w) of the irreducible prop- 
agator (self-energy) S(w) is the sum of the irreducible 
diagrams that take the system from the state N to the 
state N' and 1 = 8n N' is the unit matrix with the same 
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dimension as E(w). Note that the frequency dependence 
comes directly from the Laplace transform over the time 
r, which introduces an auxiliary line with energy hui. All 
diagrams located between the times t = r and t = 
therefore depend on u>. 



The diagrams entering the self-energy to this order are 
all drawn in Fig.|a)-h). 



B. Vertex Corrections 



A. The self-energy E(w) 

The Dyson equation allows us to appropriately sum 
up the diagrams to a certain order by calculating the self 
energy to that order and then inserting it into Eq. (gj). 

As the reservoirs are assumed to be in local equilib- 
rium, we chose to include only diagrams containing at 
most one electron-hole excitation at any given time. This 
approximation corresponds to keeping the irreducible di- 
agrams where any vertical line cuts at the most one in- 
ternal line. This is equivalent to the sequential tunneling 
approximation leading to the Master equation of ortho- 
dox SET-theory. 



To calculate the noise spectral density we also need to 
sum all vertex corrections to the same order. 

As the sequential tunneling approximation only in- 
cludes terms with at most one tunneling event at a time, 
all the vertex corrections that enter are those drawn in 
Fig. @. 



Main Result 



All diagrams which enter Eq. (24) within our approx- 
imations are drawn in Fig. Adding them together 
gives 



Sv(co) = %J2 p " Re \ NAN - 1 [( N + 1 ^U") + Nri(-"))](^N + i-A N y 



C 2 



N 



I 

U! 



[(N - l) 7iV H + Nyjfi-u)] (A N _t - A N ) 



(35) 



where = J2 N , N'Uw^fa) and P$ is the steady- 
state probability of there being N extra electrons on the 
island. In Eq. (|35| ) we have used that the steady-state 

t£(°) ust 



probabilities fulfill the relation P^ +1 



ps 



(see 



for instance 

The first terrpJn Eq. ( |35| ) corresponds to the result in 
Eq. (27) of Ref £3 but with frequency dependent tunnel- 
ing rates, while the other terms originate from the vertex 
corrections. The asymmetry in the noise with respect to 
positive and negative frequencies arises solely from these 
vertex corrections. In the limit ui — > Eq. itiq) coincides 
with the zero frequency limit given in ReftST Note that 
the frequency independent terms in the vertex corrections 
cancel. 

Using the relation Re{IlM} = -iS(u;)Im{n(w)} 
(this is the real part of Eq. (f33|) , using that the free prop- 
agator and the first order self energies are purely imag- 
inary) Eq. (|§5|) can be rewritten in matrix form, which 
gives the main result in this paper, 



E(w) 



I'M 



i* S \ (36) 



where j(u>) is a tridiagonal matrix whose elements 
are 7jv,jvM = 7wM - 7jvM, 1n+i,n(u) = 7jvM and 



7JV-i,Jv(w) = -7jvH- 

Note that contrary to E(w), j(u>) is not symmetric in 
frequency and for negative frequencies larger than the 
maximally available energy from the SET j(uj) is analyt- 
ically zero while for positive frequencies, j(uj) does not 
tend to zero but to a finite value. 



V. BACK-ACTION DURING MEASUREMENT 



When measuring with the SET, the bias voltage is typi- 
cally large enough to allow for a DC-current through the 
SET but not much larger, which implies that only the 
charge states and 1 have a non-zero steady-state prob- 
ability. In this case they are given by Pq* — JxJfi) 

and Pf = 



7r(0)+7 + (0) 

. We assume < 0, > and 

7i (0)+7o (°) 

|Aq I > |A^| so that electrons typically tunnel from the 
left to the right, and there is a finite DC-current through 
the SET. We will use these assumptions about the bias 
throughout the remaining part of the paper, except in 
the section about the off-state noise of the normal state 
SET. 
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FIG. 10: The sum of all diagrams within the sequential tunneling approximation. The propagators Iljvjv', which is the sum 
over all the diagrams that take the system from the state N to the state N' , are calculated using a Dyson equation (drawn 
graphically in Fig. H). 



A. Low-frequency regime 

In the low-frequency regime, defined as the regime 
where Jq(uj) and 7+ (ui) are exponentially small, the 
charge states and 1 are the only states energetically ac- 
cessible, also taking into account the externally available 
energy hu. In this case the matrix inversion in Eq. j3^ ) 
is easy to calculate analytically and the noise spectral 
density is given by 



S v (uj) 



2e 2 



p st 7o + H+pf 7 r(^) 



C 2 + [ 7 +H + 7o + (-^) + 7fH + 7f(-^)] 2 ' 

This expression has a very simple form: The sum of the 
steady state probabilities weighted by the inelastic tun- 
neling rates for transitions away from the state, normal- 
ized by a denominator containing the finite lifetimes of 
the states. For zero frequency this corresponds to clas- 
sical telegraph-noise. Note that Eq. (p^) is valid both 
in the normal and superconducting states, the difference 
only entering in the expressions for the rates 7^ (u>). 



B. High-frequency regime 



bias to be small enough to keep only the steady-state 
probabilities Pq* and Pf* non-zero, Eq. (|3^) gives 



Sy(w) 



2e 2 P st [ 7o + (w)+7iTM] + Pf hi M + 7i + (w)] 



C 2 



(39) 

In the high frequency limit, Tiuj 3> {Ec, eV}, all rates 
are similar and they are proportional both to the normal 
state tunnel conductance and the frequency 



±( \ ~ ±r \ huj 
7o M ~ 7i H = ^2 



1 

P£ 



1 



O(l) 



(40) 



where O(l) indicates a bias-dependent constant. This is 
valid both in the normal and superconducting states. It is 
clear that inserting Eq. @) into Eq. (|9|) gives Eq. (|||). 
It might be interesting to note that it is enough to in- 
clude four charge states to recover the full Nyquist noise. 
If only two charge states were included, an extra charge 
on the island would prevent further electrons to tunnel 
until the extra electron has left the island, and the cor- 
relation effectively would reduce the noise to that of a 
single junction. For similar reasons, for intermediate fre- 
quencies the noise should be reduced, compared to the 
Nyquist noise. 



In the high frequency limit the spectral noise density 
of the SET should be independent of the bias and be 
dominated by Jthe Nyquist noise, which in this regime 
(hu > k B T) isED 



2huR 



(38) 



X+(uR\\C) 

where Z(ui) is the impedance of the SET island to 
ground and R\\ = (l/P£ + 1/Pr) ~ • In this limit 
<C w and the matrix inversion in Eq. ( |36| ) can 
be Taylor expanded and approximated by the first term 

Still assuming the voltage 



C. Normal State SET 

For an SET operated in the normal state, the den- 
sity of states can be assumed to be energy independent 
when calculating the tunneling rates in Eq. (|2q). Using 
Pi,n(E) = p rtn (E) = p j\[ , the tunneling rates t^cj) can 
be written 



(A r N - hu) 



7atH 



(A 



N-l 



Tiuj) 



(41) 



(42) 



10 



SET emits 



SET absorbs 



-K\ -K\ o m\ K\ 



max^J 



hco 



FIG. 11: A schematic picture of the different processes in the 
SET for different frequencies. Note that negative frequencies 
correspond to processes where the SET emits energy, while 
positive frequencies correspond to processes where the SET 
absorbs energy. 



where a T n 



\ 2 P% 



is the di- 



mensionless conductivity (Rk is the quantum 
resistance and R r T is the tunneling resistance 
of junction r), T~(E) = E/(l - exp(-0E)), 
T+(E) = Eexp(-/3E)/(1 - exp(-0E)) and (3 = l/k B T. 
Note that includes both charging and biasing 

energies. 

At zero temperature, the r become step func- 
tions multiplied by a linear term, T~(E) = \E\0(E), 
T + (E) = \E\9(—E). In this limit the rates are easy to 
analyze. 

For frequencies of small magnitudes \hu\ < IAH-J 
only the charge states |0) and |1) are energetically al- 
lowed and we can use Eq. ( p7| ) to calculate the noise 
spectral density. Even though the expression in Eq. J37] ) 
looks very much like the classical expression with fre- 
quency dependent rates, this frequency dependence of 
the rates changes the behaviour quite drastically. The 
spectral noise density is no longer symmetric with respect 
to to, and there is a finite maximum energy available for 
emission from the SET, which can be seen in Fig. [l2] as 
Sy — for large negative frequencies hu; < — |Aq |. This 
means that if the energy splitting of the qubit is larger 
than the energy gained by putting an extra electron on 
the island (Aq ), there isn't enough energy available from 
the SET to excite the qubit, and the SET behaves as a 
passive load, only able to absorb energy. 

The preference of the SET to absorb energy rather 
than emit is also clear as SV( W ) > Sy(— w) for any lu > 0. 
This means that any two-level system with finite energy 
splitting driven to steady-state solely by the SET will not 
have an equal steady-state probability of both states. 



5 0.40 
3 0.20 



E/5 
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FIG. 12: Spectral noise density for an SET run in nor- 
mal mode. The calculation was done for zero temperature 
and with symmetric tunnel junctions Rr = Rl = 21.5 kfi. 
The DC-current through the SET was 1.5nA, n x = 0.25, 
E c = 2.5 K. 
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FIG. 13: The Coulomb staircase of an SCB driven to steady 
state by the SET run in normal state. The parameters we 
have used for the SCB are E qb = 2.5 K, Ej = 0.1 K, 
A = 2.5 K and for the SET we have used E c = 2.5 K, 
Rr = Rl= 21.5 kQ and n x = 0.25. 



Low-frequency regime 



For low frequencies |7k<;| < |A^| when no backward 
tunneling processes are allowed, the noise spectral den- 
sity can be written 



SV(w) 



e 2 21 /e + 2ttu [P^ofi + Pfa§] 



C 2 



4[ 7o + (0) 



■7f(0)] 
the 



(43) 



where the first term in 

/ = 2 e7o +(0) 7 r(0)/[7o + (0) + 7f(0)] 
though the SET. We see that the difference compared 



numerator 
the DC-current 



with classical telegraph-noise is the term linear in lu in 
the nominator. This is a quantum mechanical correction 
originating from the vertex corrections. In this regime 
the frequency dependent part of the tunneling rates in 
the denominator cancel. 

In the symmetrised noise Sy m (uj) = Sv(u) + Sv{— w), 
the linear term in the numerator cancels out. Thus in this 
region, for quantities that are proportional to the sym- 
metrized noise, such as the mixing-time (see Eq. (|14|)), 
the classical telegraph-noise give the same result. But for 
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other quantities, such as the steady-state probabilities of 
a qubit driven by the SET (see Eq. © and Eq. (EH) 
the difference is evident even for small frequencies. 



2. Coulomb staircase 

Using the SET to measure the average charge of the 
Cooper-pair box qubit it is reasonable to assume that the 
back-action from the SET is the dominant noise source. 
At the degeneracy point of the qubit, the energy splitting 
between its two eigenstates is Ej. If Ej < Aq we can 
use Eq. ( filf ) to calculate the Coulomb staircase (Eq. |l2] ) 
close to the degeneracy as 



(Q) = e 



8E c Tra 
1 + —ry-, — on g 

hl/e 



AE C 



, (44) 



where 5n g is the deviation from the degeneracy point 
(n g = 1/2) and we have assumed symmetric junctions 
(oq = otQ = ao or — R^p = Rt) and a symmetric 
voltage bias in the SET. Thus, close to the degeneracy, 
we will always get a linear charge increase for suitable 
choice of SET bias. In this regime the derivative is thus 
determined by the current through the SET rather than 
the Josephson energy in the qubit. 

Away from the degeneracy point, when the energy 
splitting of the qubit is increased, the low-frequency re- 
quirement for Eq. (|43| ) may not be fulfilled. In order to 
calculate the influence from the noise on the Coulomb 
staircase for arbitrary qubit gate voltage we have to in- 
clude the full expression from Eq. (35). The result for 
a typical setup is plotted in Fig. [13[ demonstrating that 
the back-action noise from the SET introduces additional 
smearing of the Coulomb staircase. 

This can be compared to the results by NazarovS, 
where the influence of the back-action of an SET in the 
normal state is calculated on a small metallic island in 
the normal state. 



3. Mixing time 

Using the tunneling rates in Eq. j4]]j4^ ) and insert- 
ing them into Eq. ( |35| ) we can calculate the mixing time 
T m ix due to the voltage fluctuation on the SET-island as 
a function oi-jthe energy splitting. Using a state-of-the- 
art RF-SETlia coupled to a qubit with realistic parame- 
ters (see caption), as shown in Fig. |lj, this would give 
a mixing time of approximately 10 fj,s. This should be 
compared with the measurement time t ms needed to re- 
solve the two charge states in the same setup which is 
about 0.4/xs. The resulting signal-to-noise ratio (SNR) is 
SNR = \J T mix / t ms 5, which indicates that single-shot 
read-out is possible. 




FIG. 14: The mixing time of an SCB caused by the SET. The 
inset shows an expanded view around Tiw ~ E q b. Thus, for an 
energy splitting of the qubit of approximately E q b, the mixing 
time is around 10 fj,s. The parameters used were Ec = 2.5 K, 
E qb = 0.8 K, Ej = 0.15 K, R R = R L = 21.5 kQ, k = 0.01, 
Idc = 9.6 nA. 



4- Off-state noise - Qubit reset 

One property of the SET used as a charge qubit read- 
out device is that it may be switched off by lowering the 
driving bias so that sequential tunneling is no longer pos- 
sible, i.e. both < Ag and < Aq . In this regime the 
voltage noise is determined by co-tunneling processes^. 
Since co-tunneling is a second order process in the tiiifc. 
neling conductance the voltage noise in the off-stateclEJ 
is several orders of magnitude smaller than the on-state 
noise. 

Taking energy exchange with the qubit into account 
there may be a first order tunneling event in the SET, 
even though the driving bias is too small for sequential 
tunneling. The energy taken from relaxing the qubit may 
stimulate a photon-assisted first-order tunnel event in 
the SET. At zero temperature the condition for such an 
event is simply AE > min{AQ, A^}. The voltage noise 
spectral density of the SET in the off-state is shown in 
Fig. [l5|. The curve has been calculated using Eq. j3^), 
with P,f = 1. 

This implies that in order to benefit from the low 
voltage fluctuations in the off-state the SET should be 
switched off by switching both the driving bias to zero 
and using the SET gate voltage to put it far into the 
Coulomb co-tunneling regime, i.e n x w 0. 

The nonlinearity of the voltage noise spectral density 
may also be used for fast relaxation of the qubit, i.e. as 
a qubit reset button. If the gate voltage of the SET is 
such that AE > |Aq| ps A^, and the driving bias is 
zero, the qubit relaxation rate is first-order in the tunnel 
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conductance, while the excitation rate is given by co- 
tunneling. The normal state SET may thus be used for 
qubit reset, or in other words as a switchable dissipative 
environment to the qubit. 
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D. Superconducting SET 



Compared with a normal state SET (NSET) the su- 
perconducting SET (SSET) shows two main differences. 
The density of states in the reservoirs is changed by the 
superconducting energy gap A, and in addition to quasi- 
particle tunneling also Cooper pair tunneling may occur. 
We will consider an SSET biased so that sequential quasi- 
particle tunneling is allowed, and in this regime Cooper 
pair tunneling may be neglected. Thus the same model as 
before can be used, only taking into account the changed 
quasiparticle density of states. As we are interested in 
an SET made out of aluminum, we use the BCS density 
of states 



p{E) = p N 



\E\ 



E 2 - A 2 



|£|-A), 



(45) 



FIG. 15: Noise spectral density of an SET in the off-state. 
Note that only contributions from positive frequencies remain, 
as no energy can be emitted from the SET within the sequen- 
tial tunneling approximation. The noiseless region is given by 

An - TlUJ > 0. 



where pjv is the density of states of the normal state. 
Inserting these into the expression for the tunneling rates 
in Eq. ( p6[) we get for zero temperature (see e.g. RefE3) 



7T a r 9{^Uuj ± eV r 

h 



2A) 



2A =p ± eV r 



-(2AThtu±eV r ) 2 \ K( 



(Hw - eV r ) 2 K( 



hu> - eV r ± 2 A , 
hui - eV r =F 2A ' 



, hoj - eV r ± 2 A 
1 Tilu - eV r i 2A 



)-E( 



hu) - eV r ± 2A 
hu - eV r T 2A 



where K(x) and E(x) are elliptic integrals of the first and 
second kind. These rates behave just as the IV curve for 
an SIS-junction. The singularities in the superconduct- 
ing density of states introduce discontinuities into the 
tunneling rates. These discontinuities will also introduce 
discontinuities in the noise spectral density. 



1. Comparison between an SSET and an NSET 

Comparing the noise spectral density of an NSET and 
an SSET (see Fig. |l6|) is not completely straightforward 
as the SSET requires considerably higher voltage bias in 
order to get sequential quasiparticle tunneling through 
the SET, i.e. \eV L -eV R \ > 4A + £' C (1 -2n x ). Therefore 
when comparing these two in the on-state (i.e. while mea- 
suring), we use the same tunnel conductance and gate 
voltage, and then choose a voltage bias that gives the 
same DC-current through the two SETs. This is moti- 
vated by the fact that the zero frequency noise is deter- 



mined by the DC-current through the SET, this biasing 
therefore yields the same zero frequency telegraph noise 
for both the SSET and the NSET. 

Apart from the discontinuities in the spectral density 
of the SSET, the finite frequency noise differs in an- 
other important aspect. Although the two SETs carry 
the same DC current, the processes producing that cur- 
rent are qualitatively different. In the superconducting 
SET biased just above the threshold the energy gain in 
each single tunnel event is quite small, determined by 
approximately max{|eVi|, |eVR|} — 2 A. The relatively 
large current is an effect of the divergent density of state 
peaks in the reservoirs. In the normal state SET carry- 
ing the same current the maximum energy that may be 
extracted from a single tunneling event is instead quite 
large, proportional to max{|eVL|, |eVn|}. 

Comparing the voltage noise spectral density for nega- 
tive frequencies, capable of exciting the measured sys- 
tem, we find that the SSET noise is zero for Huj < 
-{max{\eV L \,\eV R \} - 2A), while the NSET spectrum 
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FIG. 16: Comparing the Spectral noise density for a super- 
conducting and a normal conducting SET. The parameters 
used were R L = R R = 21.5 kfi, E c = 2.5 K, n x = 0.25 and 
Ice = 9.6 nA. 



extends down to hid w — max{|eVL|, |eVn|} « — 2A. 

Measuring the Coulomb staircase with an NSET and 
an SSET biased to the same DC-current will thus give 
different results. The Coulomb staircase is sharper for the 
SSET because the lower amount of energy extractable 
from the SET reduces the excitation rate for the two- 
level system, and the discontinuities in the noise spectral 
density of the SSET are also clearly visible, as seen in 
Fig. fL7|. Even though this is a completely different bias 
regime, similar structure appears in RefH3. 
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FIG. 17: Comparison of an ideal Coulomb staircase and a 
staircase where the qubit is driven to steady state by either 
an SET in the normal or in the superconducting state. We 
have use the same parameters as in Fig. [l3|. 



Note that the staircases in Fig. [l7| has been calculated 



for zero temperature and for a fixed voltage bias across 
the SET, and that the DC-current is different in Fig. [l7| 
and Fig. H 

When calculating the total mixing time, the sum of re- 
laxation and absorption rates enters, and the difference 
between an SSET and an NSET diminishes. The lower 
tendency for the superconducting SET to emit is com- 
pensated for by an increased tendency to absorb energy. 



03 

1 4 
'a 







- NSET 




SSET 

















4 6 
hoo/E,. 



8 



10 



FIG. 18: The mixing time due to the noise from the SET run 
either in the normal state or in the superconducting state, 
using the same DC-current through the SET (approximately 
10 nA). 



Since the mixing time due to an SSET dependends on 
the sum of the contributions from absorptive and emis- 
sive processes, it is thus not very different from an NSET 
carrying the same DC-current. An example can be seen 
in Fig. |l|. 

VI. CONCLUSIONS 

We have calculated the full frequency spectral den- 
sity of voltage fluctuations in a Single Electron Tran- 
sistor (SET), used as an electrometer biased above the 
Coulomb threshold so that the current through the SET 
is carried by sequential tunneling events. We take the 
energy exchange between the SET and the measured sys- 
tem into account using a real-time diagrammatic Keldysh 
technique. We find simple analytical expressions for the 
noise in the low- and high-frequency regimes and in be- 
tween we calculate the noise numerically. The complexity 
of the numerical calculation is limited to the inversion of 
a NXN matrix where N is the number of charge states 
involved, typically N < 5. 

Previous expressions for the voltage fluctuations, 
where the energy exchange is not taken into account, are 
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by definition symmetric with respect to positive and neg- 
ative frequencies. We show that there is an asymmetry, 
technically arising from the first order vertex corrections 
of the external vertices, so that the noise capable of ex- 
citing the measured system is always less than the noise 
that will relax the measured system, at any given fre- 
quency. The importance of this difference is shown by 
calculating the Coulomb staircase of a Cooper pair box, 
as measured by the SET. Interestingly the difference has 
a tendency to cancel in the expression for the symmetric 
noise, i.e. the sum of the positive and negative frequency 
noise. This implies that the classical calculation is a rea- 
sonably good approximation for that quantity. 

The divergence in the superconducting density of 
states results in discontinuities in the voltage noise spec- 



tral density of the superconducting SET (SSET). Com- 
pared to a normal state SET carrying the same DC cur- 
rent the SSET also has considerably less ability to excite 
the measured system. 
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